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Abstract 



Q^, We prove upper bounds on the rate, called "mixing rate", at which the von Neumann 

entropy of the expected density operator of a given ensemble of states changes under non- 
^ I local unitary evolution. For an ensemble consisting of two states, with probabilities of p and 

^ , ' 1 — p, we prove that the mixing rate is bounded above by 4a/p(1 — p) for any Hamiltonian 

of norm 1. For a general ensemble of states with probabilities distributed according to a 
random variable X and individually evolving according to any set of bounded Hamiltonians, 
■ we conjecture that the mixing rate is bounded above by a Shannon entropy of a random 

^2 ■ variable X. For this general case we prove an upper bound that is independent of the 

OO ! dimension of the Hilbert space on which states in the ensemble act. 

m 

^ ! 1 Introduction 

The problem addressed in this paper is, given an ensemble of states, S, to find an the upper 
^ ' bound on the rate, A{S), at which the von Neumann entropy of the expected density operator 
of this ensemble changes under non-local unitary evolution. The conjecture is known as 'Small 
Incremental Mixing' for an ensemble consisting of two states, S2, and it states that the mixing 
rate is bounded above by a binary entropy S{p) = —plnp — {1 — p) ln(l — p), where p and 1 — p 
are the probabilities of the two states in the ensemble. The problem, to our knowledge, was first 
introduced by Bravyi in [3]. 

We prove. Theorem \2.2\ that for any enesemble consisting of two states, the mixing rate is 
bounded above by 

A(^2) < Wvi^-P)- 
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This bound has a shape similar to that of the binary entropy, which appears in the conjecture, up 
to a factor of 2. But, unfortunately, our behavior near p = is significantly worse than pin p. 

Bravyi proved [3] the Small Incremental Mixing conjecture for a special case in which the 
expected density operator has at most two distinct eigenvalues, of arbitrary multiplicity. He 
could bound the mixing rate by 6 times the binary entropy S{p). For a general case he gave the 
dimension-independent and p-independent bound of 2, which was the best dimension-independent 
bound untill now. See Section 2 for the discussion of this conjecture. 

In our paper, see Section 3, we generalize the problem to the ensemble of any number of states, 
S, not only two, and conjecture that the upper bound should be a Shannon entropy of the random 
variable X, according to which the probabilities are distributed in the ensemble, 

A(£') < S{X). 

We prove that the mixing rate has an upper bound independent of the dimension of the Hilbert 
space the states act on. See Theorem 13.21 for the formulation of the result and Section 4 for the 
proof of the upper bound. 

Bravyi introduced the Small Incremental Mixing problem as a generalization of the 'Small 
Incremental Entanglement' conjecture, [3]. According to Bravyi, the latter conjecture was first 
proposed by Kitaev in a private communications to him. It bounds the rate of change of an 
entanglement between two parties when the system evolves under a non-local unitary evolution. 
The conjecture states that the upper bound is clnrf, where d is a dimension of a system of either 
party and c is a constant independent of either dimension. See Section 2 for the discussion of this 
conjecture. The proof of the conjecture with the constant c = 4 can be found in [6l [7]. One would 
hope to improve the constant 4 in this proof. 

The question of bounding a mixing rate by a binary entropy for an ensemble of two states is 
still open. In the special case discussed by Bravyi [3], one would hope to improve the constant 6 
in front of the binary entropy. A conjecture of bounding a mixing rate by a Shannon entropy for 
a general ensemble is open as well. 

The paper is organized as follows. In Section 2 we discuss original Small Incremental Mixing 
problem (for the ensemble consisting of two states). Small Incremental Entangling problem, the 
relation between the two and the progress on both problems. In Section 3 we generalize Small 
Incremental Mixing to a general ensemble consisting of any number of states, pose a new conjecture 
and provide our main result on the upper bound of the mixing rate. Theorem 13.21 Section 4 
contains the proof of Theorem 13.21 

Acknowledgments. We are grateful to Frank Verstraete for making us aware of the mixing 
rate problem and for his encouragement. 

2 Preliminaries 

Let T-L denote a D-dimensional Hilbert space. Let S2 = {(p, pi),(l — p, P2)} be a probabihstic 
ensemble of two states acting on "H. The expected density operator of this ensemble is a convex 
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combination p = ppi + (1 — p)p2- For any Hamiltonian H (self-adjoint operator on "H) we can 
define a time dependent state 



iHt „ JHt 



pit) =ppi + {l-p)e '"'p2e 



That is H acts locally on p2, but not on pi. 
The von Neumann entropy of this state is 



5(p(t)) = -Tr(p(t)lnp(t)). 



From the basic properties of von Neumann entropy, the following holds. 

Small Total Mixing. (Binary case) For any fixed ensemble £2, the entropy of a state p{t) 
at any time t satisfies 

SiS2)<S{p{t))<S{S2) + Sip), 

where 8(82) = pS{pi) + (1 — p)S{p2) is the average entropy of the ensemble and S{p) = —plnp — 
(1 — p) ln(l — p) is a binary entropy. 

The inequality is proved in Chapter 3 for a general ensemble of any number of states, see (13. ip . 

The analogue of the small total mixing for infinitely small times is formulated in terms of a 
mixing rate. 

A mixing rate is defined as 

dS{p{t)) 



HS2,H) 



dt 



t=o 

2.1 CONJECTURE. (Bravyi |3]) Small Incremental Mixing. 

For any ensemble £2 = {{p, pi), (1— p, P2)}, the maximum mixing rate is bounded above by a binary 
entropy. 

A{S2) : = max{|A(^, H)\ : \\H\\ = 1} 

< S{p) = —plnp — {1 — p) ln(l — p). 

Some useful formulas for the mixing rate are 

A{S2,H) = -ipTT{[pi,\np]H) (2.1) 
= i{l-p)TT{[p2,\np]H), 

and 

A(£'2) =pTr|[pi,lnp]| 
= p||[p2,lnp]||i, 

here the maximum is achieved for H = 1 — 2R, with R being a projector on the negative eigenspace 
of z[pi,lnp]. The norm || ■ ||i is a trace-norm. 

Bravyi [3] proved that ^{82) < 6S{p), where p has at most two distinct eigenvalues of arbitrary 
multiplicity. 

Our result for an ensemble of two states is the following theorem. 
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2.2 THEOREM. (Binary case) For any binary ensemble £2 = {(p, Pi), (1 — P,P2)}, the max- 
imum mixing rate is bounded above 

A{£2) < A^/p{l-p). 

The proof of this theorem is given in Chapter 3 for a more general case, when the ensemble 
consists of any number of states, see Theorem 13.21 

Although we do not pursue this direction, we mention the following question posed by Aude- 
naert and Kittaneh pQ. 

2.3 PROBLEM. (Audenaert, Kittaneh |T]) Let A and B be arbitrary positive semi-definite 
d X d matrices with a = TrA and b = TrB. For what functions / : M — >■ M does there exist a 
constant c, independent of d, A and B, such that 

II [B, f{A + B)] 111 < c(F(a + 6) - F(a) - F{b)), 

where F{x) = f{y)dy? 

Note that if A = ppi, B = {1 — p)p2 and /(x) = Inx, Problem 12.31 becomes Conjecture 12. 1[ 
As mentioned in the introduction, the Small Incremental Mixing problem is a generalization of 
the Small Incremental Entangling conjecture. To formulate the later conjecture, we suppose that 
two parties, say Alice and Bob, have control over systems A and B. Both systems evolve according 
to a non-local Hamiltonian Hab- In time entanglement between A and B can be generated. In 
ancilla-assisted entangling both parties have access to additional subsystems, called local ancillas, 
i.e. Alice is in control of two systems A and a and Bob is in control of B and b. Alice and Bob 
start with a pure state p(0) = |^) (\1/| on the system aABb. 
A time dependent joint state of Alice and Bob is 

p{t) = U*it)\^) {^\Uit), 

where U{t) = (8> e^^^^t /^^ is a unitary transformation. The joint state of Alice and Bob is pure 
at any time. 

One of the ways to describe the entanglement between Alice and Bob is to calculate the 
entanglement entropy 

E{p{t)) := S{paA{t)) = -TTPaA{t)\npaA{t), 

where paA{t) = Ti^BbP{t) is a state that Alice has after time t. Since the joint state is pure, the 
entanglement entropy is also can be calculated from the state that Bob has E{p{t)) = S{pBb{t))- 

Small Total Entangling. The total change of the entanglement E{p(t)) is at most 2lnd, 
where d = min{dim(y4), dim(i?)}. See [2] for the proof. 

A problem of bounding the infinitesimal change of the entanglement is formulated using the 
entangling rate. 

The entangling rate is defined by 

dE(p(t) 



t=0 
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After calculating the derivative, the entangling rate can be expressed as 

r(^, H) = -iTlfHABlPaABMiPaA) ® Ib] 



= -lTr(^HAB[PaABM{PaA ® ^■^^^^ )]). (2.2) 

Similarly, r(^, H) = -iTi(^HAB[pABb, Ia ® HpBb)]) 

2.4 CONJECTURE. (Bravyi P) Small Incremental Entangling. 

There is a universal constant c such that for all dimensions of ancillas a, b and for all states |\&), 
the following holds 

r(*,//) < c\\H\\\nd, 

where d = min{dim(y4), dim(S)}. 

This problem was studied by many authors. For the case when A and B are quibits, Childs et 
al |1] give upper bounds for the entangling rate and show that they are independent of the ancillas 
a and b. Wang and Sanders [8] proved that T{H) < /3, where /3 ~ 1.9123, for an uncorrelated 
Hamiltonian H = Ha ® Hb, when Ha{b) = H'^^^y Child et al [5] also proved an upper bound 
T{H) := maxij, r(\E', if) = [3\^a^b ^ /^ll-^ll for the ancilla- assisted case and for an arbitrary 
uncorrelated bipartite Hamiltonian H = Ha ® Hb, where {^b) is the difference between the 
largest and smallest eigenvalues of Ha (Hb). For an arbitrary bipartite Hamiltonian Bennet et al 
[2] proved that T{H) < cd'^WHW, where c does not depend on a or b. Finally, the conjecture was 
proved by Hutter and Wehner [6], [7], where they showed that for any initial state p(0) (pure or 
mixed) and any Hamiltonian H, r{p{0),H) < 41n(i||if||. 

Bravyi [3] proved that Small Incremental Mixing with a constant c in front of the Shannon 
entropy implies Small Incremental Entangling with a constant 4c, by choosing particular ensemble 
of states: £2 = {((1 — dim{B)~'^) , HaAB) , dim{B)~^ , paAB)} ■ Here without loss of generality it was 
assumed that B < A and paAB is a state such that the expected density operator of the ensemble 
is of the form appearing in ( 12. 2p 

PaA ® ^ = (1 - dimiBy^^fiaAB + dim(5)^VaAB. (2.3) 

In Lemma 1 |3] it was proved that such a state PaAB exists. Applying (12. ip to this ensemble, 
we have that A{S2,H) = dim(i?)~^r(\l/, i/), which shows that Small Incremental Mixing implies 
Small Incremental Entangling. 

Using Bravyi's proof, our bound oi A^/p(Y^^J)j ioi the Small Incremental Mixing problem leads 
to a bound of 4(i||if|| in Small Incremental Entangling. 

3 Small Incremental Mixing for an ensemble consisting of 
any number of states 



We generalize the Small Incremental Mixing problem to an ensemble of any number of states in 
the following way. 
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Let X be a random variable with probability density px{x), i.e. the probability that the 
realization x occurs is px{x), where the realization x belongs to a set X. Consider a probabilistic 
ensemble of states £ = {px{x), px}x£X, i-e. Px is a density matrix on "H which occurs with 
probability pxix), where ^^Px{x) = 1. 

The expected density operator of the ensemble S is convex combination of density matrices 
p = J2xP^i^)Px- For collection of Hamiltonians H = {Hx}xex define a time-dependent state 



-iHxt „ „iHxt 

UxyX)ti 

xex 

Note that one of the states could always be left invariant, i.e. one of the Hamiltonians could 
always be taken as an identity /, but to simplify the notation we shall write a time evolution for 
all states. 

The von Neumann entropy of this state p{t) is 

5(p(t)) = -Tr(p(t)lnp(t)), 

which satisfies the following property. 

Small Total Mixing (General case). For any fixed ensemble £, the entropy of a state p(t) 
at any time t satisfies 

S{S)<S{p{t))<S{S) + S{X), (3.1) 

where S{S) = '^xPx{x)S{px) is the average entropy of an ensemble S andS(X) = — '^^px{x)\npx{x) 
is a Shannon entropy of a classical random variable X. 

The lower bound follows from concavity property of the von Neumann entropy and the invari- 
ance of the entropy under unitary transformation 

m = = E^^(^)^(^~^''^Vx.e^^^*) < 5(j]px(x)e-^^^Vxe^^^*) = 5(p(t)). 

X X xdX 

To see the upper bound, form a classical-quantum state 

p''^(0 = X]px(a;) \x) (x|^®p,(t), 

X 

that acts on a tensor product of a classical space X and a quantum system A represented by a 
Hilbert space %. Here, Pxif) = e~*^"=*pxe*^="*. The entropy of the classical-quantum state is 

Sip''^it)) = SiX) + SiS). 

For a classical-quantum state p^^ the relative entropy ^(XIA) := S{p^^) — S{p^) > is always 
non-negative, therefore S{p{t)) < S{p'^^(t)). This proves the Small Total Mixing property. 

The analogue of the small total mixing for infinitely small times is formulated in terms of a 
mixing rate. 

A mixing rate is defined similarly to the binary 

dS{p{t)) 



dt 



(3.2) 

£=0 
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3.1 CONJECTURE. Small Incremental Mixing. 

For any ensemble S = {{px{x), Px)}xex, the maximum mixing rate is bounded above by a Shannon 
entropy. 

A{S) : = max{|A(^,H)| : -I < < I , x e X} (3.3) 
< S{X) = - ^^pxix) Inpxix). 

X 

A mixing rate can be written explicitly by calculating a derivative of the entropy ^{p{t)) = 
-Tr(^^lnp(t)), att = 

A(^,H) = -iY,Px{x)Tt{[H^,Px] Inp) 

X 

= -tJ2Px{^)MHx[PxMp]). (3.4) 

X 

Note that in the definition of the mixing rate (13. 3p the maximum is taken over all Hamiltonians 
such that —J < < I. 

For any Hermitian operator A with Tt{A) = 0, 

max{Tr(i7yl) : -I < H < 1} = 2max{Tr{HA) ■.0< H <I}. 

This property can be easily observed by expressing H = 2R — I, where < R < I. 

Therefore the maximum in (13. 3 p can be taken over the non-negative Hamiltonians bounded 
above by identity operator. 

A(^) = 2max{|A(£,H)| ■.0<H^<I,xe X}. (3.5) 

Note that, similarly to the binary case, the maximum is achieved for the Hamiltonians being 
a projector onto a positive eigenspace of z[pa;,lnp]. 

Bravyi's proof of the Small Incremental Mixing problem for an ensemble of two states and for 
a state p with binary spectrum can be easily generalized to the case of an ensemble containing 
any number of states with p still restricted to having a binary spectrum. Also the constant 6 in 
front of the Shannon entropy S{X) remains. 

In Section 3 we prove the following theorem, showing that the maximum mixing rate is bounded 
above by a constant independent of the dimension D of the Hilbert space "H that states act on. 

3.2 THEOREM. (General case) For a fixed ensemble S = {px{x), Px}x£X the maximum 
mixing rate 1^3. 3\) is bounded above 

A{£) <2J2Y1 ^Px{x)px{y) = 2 (J2 VM^) ' - 2. 

X y^x X 

For a binary ensemble S = {(p, pi),((l — p),p2)} Theorem 13.21 gives the upper bound of 
4y^p(r^^p), as claimed in Theorem 12.21 
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4 Upper bound on the mixing rate 

In this section we prove Theorem 13.21 

In eq. (13 ■4p , express the logarithm log p by the formula 

\nx= f (-^ ^—\dt. 



.l + t x + t^ 

We may assume that < < I for every x G A", as noted before f l3.5p . to calculate the mixing 
rate. Then, for reasons explained below, 

A(^,H) = -iY,Pxix)TTiH,[p,Mp]) 

1 N 

dt 



-iy^Pxix) / Tt(H^[p^,— 
Jo ^ P + 

/"°° / I 1 \ 

iJ^Pxix) / Tr -[p,,p + t]— -Wt 

V -^0 \ P + t p + tJ 

/ 1 1 \ 

Tr i^H^ [pa^ , py] j dt 

^'y.^Vx{x)Vx{v) / Tr([p^.,pj,]— -if^— -jdt. 



Here in the third equality we used that [A, B~^] = B~^[A, B]B^^. In the fourth equality we wrote 
p as a convex combination of states py and eliminated the commuting terms. In the last equality 
the cyclicity of the trace is used. 

For any < < /, we have < (p + t)~^H{p + t)'^ <{p + t)~^. 

Therefore, continuing our calculations, with := f^{p + t)~^Hx{p + t)~^dt < p~^, x E X 

I / 1"°° 1 1 

|A(^,H)| = 5^5^px(xK(?/)Tr [p,.,p,] / ^M—-dt 

X yjLx Jo P + T^ P + r 

< ^^Px{x)px{y)(\^^pxpyK:,\ + iTrpyp^Ji, 

X y^x 

^^^Px {x)px {y) (Tr I \/lC,pxpy \fK'x \ + Tr | \fK':^pyp^ \fK, 

X y^x 

< 2 E E (^)P^ ^J'^<PlKx)Tr{plK,) 



X y^x 



< '2^^ VPx{x)px{y)\ ^r(p^{px{x)p^)p ^)TT(py{px{y)py)p 1 



X y^x 
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In the first inequality we put the absolute value inside the sums, wrote the commutator [p^^ py\ = 
PxPy — PyPx and used the triangular property of the absolute value. The second inequality follows 
from the cyclicity of the trace and by moving the absolute value inside the trace. The third 
inequality follows from Cauchy-Schwartz inequality for traces: 



The fourth inequality follows from the upper bound on < The fifth inequality is obtained 
from the definition of p as a convex combination of non-negative density operators p^, therefore 
p > Px{x)px for any x E X, and Trp^; = 1, x E X . 

From (13. 5p we obtain an upper bound for the mixing rate 
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Tt\AB\ < ^/Tt{A*A)Tt{B*B). 
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